Fast collisionless reconnection and electron heating in strongly magnetized plasmas 
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Magnetic reconnection in strongly magnetized (low-beta), weakly collisional plasmas is investi- 
gated using a novel fluid-kinetic model [Zocco & Schekochihin, Phys. Plasmas 18, 102309 (2011)] 
which retains non-isothermal electron kinetics. It is shown that electron heating via Landau damp- 
ing (linear phase mixing) is the dominant dissipation mechanism. In time, electron heating occurs 
after the peak of the reconnection rate; in space, it is concentrated along the separatrices of the 
magnetic island. For sufficiently large systems, the peak reconnection rate is cE m ax ~ 0.2vaB Vj o, 
where va is the Alfven speed based on the reconnecting field B y ,o. The island saturation width is 
the same as in MHD models except for small systems, when it becomes comparable to the kinetic 
scales. 

PACS numbers: 52.35.Vd, 96.60.lv, 52.35.Py, 52.30.Gz 



Introduction. Magnetic reconnection is a reconfigura- 
tion of the magnetic field in a plasma via localized un- 
freezing of the magnetic flux [1]. It is commonly associ- 
ated with energy release in astrophysical and laboratory 
plasmas, but the details of the energy conversion mech- 
anisms and partition between particle species and fields 
are poorly understood. This Letter focuses on a key as- 
pect of this issue: the conversion of the magnetic energy 
into electron internal energy, i.e., electron heating. 

We consider a fundamental reconnection paradigm: 
the tearing mode in a periodic box [2]. The tearing in- 
stability leads to the opening, growth and saturation of 
a magnetic island [2-6]. Since the saturated width of 
the island is in general macroscopic, it cannot depend on 
the microphysics of the plasma and, in particular, should 
not depend on its collisionality (as we will show). This 
implies that the total fraction of the initial magnetic en- 
ergy converted to other forms of energy from the be- 
gining to the end of the evolution of the tearing mode 
must be the same in collisional and collisionless plasmas. 
In a periodic (closed) system, energy cannot be lost via 
bulk plasma outflows. Therefore, the magnetic energy 
difference between the initial and final states must be 
accounted for by conversion into the thermal energy of 
the particles. In collisional plasmas, this is achieved by 
Ohmic and viscous heating [7]. However, many natu- 
ral systems where reconnection occurs are only weakly 
collisional; the only available heating channels then are 
Landau damping and electron viscosity, both of which 
ultimately rely on the electron collision frequency being 
finite, though possibly arbitrarily small. In this Letter, 
we show that electron heating via linear phase mixing as- 
sociated with Landau damping is the main energy conver- 
sion channel in weakly collisional reconnection in strongly 
magnetized (low-beta) plasmas. 

Equations. We use a fluid-kinetic approximation ap- 
plicable to low-beta plasmas ("KREHM" [8]). Its main 



distinctive feature, for our purposes, is the coupling of 
Ohm's law to the electron (drift) kinetic equation via the 
electron temperature fluctuations. The kinetic equation 
allows for both collisions and Landau resonance, thus en- 
abling "collisionless" electron heating. 

We work in the low-/3 regime, where /?, the ratio of 
the plasma to the magnetic pressure, is ordered similar 
to the electron- ion mass ratio m e /rrii. Let us define the 
perturbed electron distribution function to lowest order 
^/m e /rrii ~ and in the gyrokinetic expansion [9- 
9e + (Sn e /noe + 2v\\U\\ e /v^ he )F 0e , where 
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12] as Sf ( 

Fo e is the equilibrium Maxwellian, vth 



y / 2T 0e /m e 

is the electron thermal speed (with Tb e the mean elec- 
tron temperature), v\\ is the parallel velocity coordinate, 
5n e /noe is the electron density perturbation (the zeroth 
moment of Sf e ) normalized to its background value no e , 
and u\\ e = (e/cm e )deVj_A|| is the parallel electron flow 
(the first moment of 5f e ] A\\ is the parallel component 
of the vector potential and d e = c/uj pe is the electron 
skin depth). All moments of Sf e higher than Sn e and 
u\\ e are contained in g e , e.g., the electron temperature 
perturbation is ST\\ e /T 0e = (l/n 0e ) / d 3 v (2vj/v^ he ) g e . 
The dynamics of the plasma is described by the evolution 
equations for 5n e /no e , A\\ and g ei which, for the 2D case 
considered here, read [8]: 
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where C[g e ] is the collision operator [13], {...,...} is the 
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Poisson bracket and d/dt = d/dt + c/ B z {(p, ...}, with B z 
the out-of-plane magnetic guide-field and cp the electro- 
static potential. The latter is obtained via the gyroki- 
netic Poisson's law [10], 5n e /no e = 1/t(Tq — l)eip/To e , 
where r = Xbi/Xb e and To is the real-space operator that 
is the inverse Fourier transform of To (a) = Io(a)e~ a , 
with Iq the modified Bessel function and a = k\pl/2 
{pi = Vthi/^i is the ion Larmor radius). 

Eq. (3) shows that the popular isothermal closure [14], 
g e = 0, is not a solution of that equation unless 
{Am, Vj_Am} = 0, a condition that cannot describe a re- 
connecting system [though it does describe the (macro- 
scopic) island saturation [4, 5], as we will find]. This 
means that at least the possibility of electron heating in 
weakly collisional reconnection cannot be ignored. 

Numerical Details. To simplify the solution of 
Eq. (3), note that it does not contain an explicit 
dependence on the perpendicular velocity coordi- 
nate, v±. If we ignore any such a dependence that 
is introduced by the collision operator, v± can be 
integrated out, so g e = g e (x,y,v\\,t). Next, we 
introduce the Hermite expansion g e (x, V-> t, vu) = 

Y,m=2 H rn(v\\/v the )g m (x,y,t)Foe(v\\)/V2 m rn ] - 
(go = 9i = because Sn e and u\\ e have been ex- 
plicitly separated in the decomposition of 5f e adopted 
above). Eq. (3) then unfolds into a series of coupled, 
fluid-like equations for each of the coefficients g m : 



dt 
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where we have adopted a model (hyper) collision opera- 
tor with Vcoii = l/(AtM 4 ), where M is the index of the 
highest Hermite polynomial kept in a simulation, and At 
the timestep. Thus, in our simulations, M is a proxy for 
the collision frequency, with higher values of M corre- 
sponding to less collisional systems (at the large values 
of M reported here, At - M -1 / 2 so v coU ~ M~ 7 / 2 ). 
Note that the more familiar (and more physical) in such 
ID problems Lenard-Bernstein collision operator would 
instead be mu e ig m [8] (y e i is the electron-ion collision fre- 
quency); because of the linear dependence on m of this 
operator an unfeasibly large number of Hermite poly- 
nomials would need to be kept in order to resolve the 
velocity-space cutoff as v e i —> 0. This is why we use 
hyper-collisions instead. 

Eqs. (1,2,4) are solved numerically using a pseudo- 
spectral code [15]. The spatial grid size is 384 2 . The res- 
olution in velocity space is set by M and ranges from 30 
to 500. Hyper-diffusive terms of the form ^Vj^, where 
vh = 0.25/At(Ax/7r) 6 , with Ax the grid spacing, are 
added to the RHS of all equations (including Eq. (4)). 
These are required to prevent the unbounded thinning 
of the current layer in the nonlinear regime [16-18], 
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FIG. 1. Time traces of: (a) the reconnection rate for different 
values of collisionality (represented by M; larger M means less 
collisions); (b) the rates of dissipation via Landau damping 
and hyper-dissipation, for the least collisional case (M = 500); 
(c) the rate of dissipation via Landau damping for different 
values of M. All runs had A a = 20. 



but will not, as we will discover, dissipate much energy. 
Physically, they stand in for collisional and collisionless 
electron finite Larmor radius effects (which are formally 
small in the KREHM approximation). 

The hyper- diffusive and hyper-collisional terms are the 
only dissipative terms we employ. We shall distinguish 
their effects: "hyper-viscous" and "hyper-Ohmic" refer, 
respectively, to the dissipation terms in Eq. (1) and in 
Eq. (2); "Landau" refers to the dissipation arising via 
both hyper-collisions and hyper-dissipation in Eq. (4). 
The rationale for the latter is that the only way that en- 
ergy can arrive at any g m with m > 2 is via phase mixing. 
Once there, how exactly it dissipates is not important: 
what we wish to investigate is the relative importance of 
dissipation via phase mixing (i.e., Landau-damping) vs. 
(hyper-) viscosity or resistivity. 

The equilibrium in-plane magnetic field is B ym = 
—dA\\ eq /dx, with A\\ eq = A\\ / cosh 2 (x/a), where a is the 
(normalizing) equilibrium scale length. The (normaliz- 
ing) Alfven time is defined as ta = cl/va with va the 
Alfven speed based on B y $ = max \B y ^ eq \ = 1. The 
simulations are performed in a doubly periodic box of 
dimensions L x x L y , with L x /a = 2tt and L y such that 
k y = 2i\ajL y yields the desired value of the tearing insta- 
bility parameter Ala = 2(5-fe 2 )(3+£; 2 )/(fe 2 y / 4 + /c 2 ) [6]. 
For simplicity and numerical convenience, we set pi = 
p s = d e = 0.25a, where p s = pi/\/2r. Thus, this study 
does not address any effects associated with the scale sep- 
aration between ions and electrons — a reasonable first 
step for low-/3 plasmas. 

Reconnection Rate. Before we discuss the electron 
heating, let us first report some of the more standard 
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FIG. 2. System configuration for our least collisional simu- 
lation (M = 500) at the time of maximum reconnection rate 
(top row) and of maximum dissipation rate (bottom row) [see 
Figs. 1(a) and (b), respectively]. From left to right, the plots 
show: A || (lines: full/dashed are positive/ negative contours) 
and j|| = — c/47rV5_A|| [colors, ranging from blue (negative) 
to red (positive)]; cp (lines) and ST\\ e (colors); the total dissi- 
pation via g m 's (i.e., via Landau damping). The horizontal 
dashed line marks the location of the domain cut where the 
distribution function of Fig. 4 (bottom row) is plotted. 



features of our simulations. The system evolves in time 
from the linear to the nonlinear stage; saturation occurs 
when all the initially available flux has been reconnected. 
The time-traces of the reconnection rate, defined as the 
value of the parallel electric field at (x, y) = (0, 0) (the 
X-point), Ex, are plotted in Fig. 1(a). As shown, the 
reconnection rate is entirely independent of collisionality 
(parameterized by M), consistent with the fact that our 
simulations are in the weakly-collisional regime, where 
the frozen-flux constraint is broken by electron inertia, 
not the collisions. The maximum value is c£?™ ax w 
0.22vAB yi o, similar to the fast reconnection rates ob- 
tained in the opposite limit of weak guide-field [19], 
and in qualitative agreement with [20, 21]. Regarding 
the dependence of £?™ ax with system size (not shown), 
we found that i?™ ax increases with A', asymptoting to 
cE™ ax ~ 0.2v A B y , for A'a > 10. 

Fig. 2 depicts the system configuration for our least 
collisional simulation (M = 500) at the time of maximum 
reconnection rate [t/rA = 22.7; see Fig. 1(a)] (top row) 
and at the time of maximum dissipation rate [t/rA = 
29.4; see Fig. 1(b)] (bottom row). A typical X-point 
geometry is seen, accompanied by a quadrupole structure 
exhibited by ST\\ e (and by 5n ei not shown) [22] [23]. 

Saturation. Fig. 3 shows the saturation amplitudes 
obtained in our simulations ("KREHM"); overplotted 
are the MHD (p Syi = d e = 0) results from Ref. [6] 
("RMHD"). Also shown is the prediction from MHD 
theory [4, 5] ("POEM"), valid in the small- A'a regime. 
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FIG. 3. Saturated flux ^ sat vs. A'. "KREHM" are the 
data from the kinetic simulations; "RMHD" are MHD results 
from [6]. The dotted line ("POEM") is the prediction from 
MHD theory [4, 5]. Inset: time-traces of flux for A'a = (1, 2). 
In both plots, the horizontal line is the saturated flux corre- 
sponding to a full island width of 2d e /a. 



As anticipated, we find that saturation in our weakly 
collisional, kinetic simulations is well described by the 
(isothermal) MHD model. At small A'a, this agreement 
breaks down because the island saturation amplitude be- 
comes comparable to the kinetic scales (d e , p s ,i)- Satura- 
tion then becomes a slow, diffusive process, as illustrated 
by the inset plot: the islands slowly expand until their 
width becomes ~ d e . Thus, there is a lower limit to 
the saturation amplitude, set by d e (cf. [24]); indeed, the 
frozen-flux condition precludes the definition of magnetic 
field lines at sub-d e scales. 

Electron heating. We have established that the 
amount of energy converted during the evolution of the 
system is independent of the collisionality; we have also 
checked (not shown) that for all but the smallest systems, 
the energy converted into electron heating is a significant 
fraction of the initial (magnetic) energy, reaching ~ 60% 
for the largest systems. We now turn to the main focus of 
this Letter: how is the energy converted. Fig. 1(b) shows 
the time traces of the dissipation rates for A'a = 20 and 
M = 500. We see that dissipation happens almost exclu- 
sively via the phase mixing/Landau channel. We reiter- 
ate that, a priori, the system is free to choose between 
different dissipation channels, i.e., the fact that phase 
mixing is the preferential dissipation mechanism is not 
hard- wired and is, therefore, a remarkable demonstration 
of the dominance of Landau damping over other dissipa- 
tion mechanisms in weakly collisional reconnection [25]. 

Another noteworthy feature in Fig. 1(b) is the time lag 
between the peaks of the reconnection and dissipation 
rates. This implies that magnetic energy is not directly 
dissipated by the reconnection process itself [26] . Instead, 
it is first converted to other forms and then dissipated. 



Indeed Landau damping dissipates the electron free en- 
ergy, / dxdy/V J dv ll T 0e gl/(2F 0e ). Fig. 1(c) shows that 
the time lag increases weakly with decreasing collisional- 
ity (see below). 

A related question is where in the domain the heating 
is occurring. The plots in the right column of Fig. 2 show 
that there is no significant heating in the current sheet; 
instead, it happens predominantly along the separatri- 
ces of the island. This is consistent with the existence 
of a time lag between reconnection and dissipation: the 
magnetic energy that is converted in the current sheet 
is mostly channeled into ion and electron kinetic energy. 
Both species flow downstream predominantly along the 
separatrices. Deceleration of these flows converts kinetic 
energy into the free energy (or entropy) of the electrons 
(plotted in the top row of Fig. 4). It is this conversion 
that constitutes the energetics of Landau damping and 
results eventually in electron heating. 

A detailed understanding of the electron heating 
process is yielded by the spectral maps [in the two- 
dimensional Fourier-Hermite (k±,m) space] of the elec- 
tron free energy and its dissipation. These are shown in 
the two top pannels of Fig. 4 for M = 500. As time 
advances, the electron free energy cascades to higher val- 
ues of m, corresponding to the formation of small scales 
in velocity space, i.e., phase mixing. During this pro- 
cess, energy dissipation occurs via the hyper-diffusive 
terms, acting at large values of k± (see middle pannel, 
second row). At later times, large enough values of m are 
reached so the collisional dissipation becomes important 
and, indeed, dominant; this ra-cutoff is clearly visible in 
the rightmost plot of the second row of Fig. 4. 

Let us estimate the velocity-space dissipation scale. 
Linearising Eq. (4) for a given k y we find that the electron 
free energy spectrum E m = \g m \ 2 /2 satisfies [8] 
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Setting dE m /dt = 2'yE rn , we obtain 
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where ra 7 = [k y B y /B z v the /(2 V^t)] 2 ? m c 
[9/(2y/2)kyB y / B z v the /is co u] 2/9 and C(k y ) is some 
/^-dependent constant. As is evidenced by Fig. 1(a), 7 
is independent of collisions and thus so is ra 7 . Therefore, 
while the mode is strongly growing, m 7 < m c and so the 
collisional cuttoff cannot be reached. This explains why 
the peak of the dissipation rate must occur later than 
that of the reconnection rate. As reconnection proceeds 
into the saturation regime, 7^0, so, regardless of 
how small v co u is, eventually m 7 > m c and, from then 
onwards, the Her mite spectrum cutoff is determined 
by m c . In our simulations, this happens at t « 26ta, 
with no appreciable dependence on collisionality (M), 
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FIG. 4. Electron free energy spectrum (top row), dissipation 
spectrum (middle row) and a cut at y/a — 2.3 (cf. Fig. 2) of 
the distribution function (bottom row) at the early nonlinear 
stage (left) and at the peaks of the reconnection rate (center) 
and dissipation rate (right) for a run with A'a — 20 and 
M = 500. 



because the decrease of 7 is fast [see Fig. 1(a)]. The 
inset in Fig. 5 shows the time lag between the peaks 
of the reconnection and dissipation rates vs. M. The 
logarithmic dependence is due to the fast decay of 7, and 
the consequent rapid increase of m 7 ~ 7" 2 to overtake 
m c , thus enabling dissipation. The weak dependence of 
the lag on collisions implies that dissipation occurs in 
finite time even for weak collisionality. 

The value of m = m pea k at which most energy is dis- 
sipated is the solution of d(v co um 4l E rri )/dm = 0, with 
E m given by Eq. (6) in the regime m 7 ^> m c . This 
yields m pea k = (9/7) 2 / 9 m c . This expression, evaluated 
for k y = 1 [27], is compared in Fig. 5 with the numerically 
determined value of m at which the dissipation peaks. 
The remarkable agreement that is obtained shows that 
the electron heating we observe is the result of linear 
phase mixing. 

The phase mixing process discussed above is illustrated 
by the plots in the bottom row of Fig. 4 showing the 
electron distribution function g e , taken at y/a = 2.3 (i.e., 
inside the island; see Fig. 2). The progressive creation 
of finer scales in velocity-space, a textbook signature of 
phase mixing, is manifest (cf. [28, 29]). 

Conclusions. This Letter presents the first detailed 
investigations of electron heating caused by magnetic 
reconnection in strongly magnetized, weakly collisional 
plasmas. Using a novel fluid-kinetic framework [8], we 
were able to show that linear phase mixing/Landau 
damping is the main mechanism for energy conversion 
and electron heating. Reconnection and electron heat- 
ing are causally related, but temporally and spatially 
disconnected: heating happens after most flux has re- 
connected, and along the island separatrices, not in 
the current sheet. Our other key conclusions are: (i) 
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FIG. 5. Value of m at which most energy is dissipated, m pea k, 
as a function of collisionality, v co ii. Inset: time lag between 
the peak of the reconnection rate (r^ ax ) and the peak of the 
dissipation rate, (r^ ax ), as a function of v co ii. 



the maximum reconnection rate in the asymptotically- 
large-guide-field limit is as fast as in the no-guide-field 
limit, cE max ~ 0.2v AByfi, provided that the system is 
large enough; (ii) the saturation amplitude in the kinetic 
(weakly collisional) regime is identical to that in MHD 
(collisional) systems [6], as long as the island is large 
compared to the kinetic scales. The electron inertia scale 
appears to provide the lower boundary on the saturation 
amplitude — this result may be important to the under- 
standing of magnetized turbulence (e.g., [30-33]), as it 
effectively sets the minimum fluctuation amplitude. 
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